c     forcing fields for momentum equations
      common /force/ uforce, vforce
      real uforce(nx,ny)
      real vforce(nx,ny)
c
     
